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CrossMark 

Abstract 

Phase  matching  of  circularly  polarized  high-order  harmonics  driven  by  counter-rotating  bi- 
chromatic  lasers  was  recently  predicted  theoretically  and  demonstrated  experimentally.  In  that 
work,  phase  matching  was  analyzed  by  assuming  that  the  total  energy,  spin  angular  momentum 
and  linear  momentum  of  the  photons  participating  in  the  process  are  conserved.  Here  we  propose 
a  new  perspective  on  phase  matching  of  circularly  polarized  high  harmonics.  We  derive  an 
extended  phase  matching  condition  by  requiring  a  new  propagation  matching  condition  between 
the  classical  vectorial  bi-chromatic  laser  pump  and  harmonics  fields.  This  allows  us  to  include 
the  influence  of  the  laser  pulse  envelopes  on  phase  matching.  We  find  that  the  helicity  dependent 
phase  matching  facilitates  generation  of  high  harmonics  beams  with  a  high  degree  of  chirality. 
Indeed,  we  present  an  experimentally  measured  chiral  spectrum  that  can  support  a  train  of 
attosecond  pulses  with  a  high  degree  of  circular  polarization.  Moreover,  while  the  degree  of 
circularity  of  the  most  intense  pulse  approaches  unity,  all  other  pulses  exhibit  reduced  circularity. 
This  feature  suggests  the  possibility  of  using  a  train  of  attosecond  pulses  as  an  isolated 
attosecond  probe  for  chiral-sensitive  experiments. 

Keywords:  circularly  polarized  high  harmonic  generation,  phase  matching,  ultrafast  chiral 
physics,  attosecond  pulses 
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1.  Introduction 


Original  content  from  this  work  may  be  used  under  the  terms 
of  the  Creative  Commons  Attribution  3.0  licence.  Any 
further  distribution  of  this  work  must  maintain  attribution  to  the  author(s)  and 
the  title  of  the  work,  journal  citation  and  DOI. 


Light  sources  based  on  high  harmonic  generation  (HHG)  are 
used  for  a  variety  of  applications  in  ultrafast  spectroscopy  and 
high-resolution  imaging  of  atoms,  molecules,  condensed 
matter  and  nano-structures  [1-14].  While  many  schemes  were 
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developed  for  controlling  the  spatial,  temporal  and  spectral 
properties  of  HHG-based  light  sources  [15-19],  until  recently 
bright  (and  useful)  HHG  radiation  was  largely  limited  to 
linear  polarization. 

The  difficulty  in  controlling  the  polarization  of  HHG 
originates  from  the  fact  that  high  harmonics  are  emitted 
through  a  coherent  recollision  phenomenon  [20-22].  An 
electron  is  first  ripped  from  an  atom  or  ion  by  a  strong  laser 
field.  The  electron  then  accelerates  forth  and  back  in  the 
oscillating  field,  and  finally  recombines  with  the  parent  ion 
while  emitting  a  high-energy  photon.  If  the  driving  laser  field 
is  linearly  polarized,  the  recollision  (and  recombination) 
probability  is  high  since  the  trajectory  of  the  electron  is  linear. 
This  ‘head  on’  recollision  yields  linearly  polarized  HHG 
bursts.  Introducing  ellipticity  into  the  driving  laser  field 
results  in  the  generation  of  elliptically  polarized  harmonics 
[23-25].  However,  in  this  case  the  HHG  efficiency  drops 
drastically  since  the  electron  trajectory  misses  the  parent  ion 
[26].  Although  many  theoretical  approaches  have  proposed 
paths  toward  generation  of  bright  highly  elliptical  or  circu¬ 
larly  polarized  high  harmonics,  the  experimental  progress  was 
quite  limited  [27-40].  A  notable  method  proposed  to  generate 
circularly  polarized  HHG  from  circularly  polarized  counter¬ 
rotating  [27,  29,  30,  41]  or  co-rotating  [42,  43]  bi-chromatic 
laser  fields.  A  single  pioneering  experimental  work  from  1995 
reported  the  generation  of  high  harmonics  using  this  geometry 
with  counter-rotating  fields  [44].  While  indication  for  the 
correct  selection  rules  was  observed  in  that  experiment,  no 
attempt  was  made  to  measure  the  polarization  state  of  the 
harmonics. 

In  2013,  Fleischer  et  al  reported  that  the  polarization  of 
high  harmonics  can  be  fully  controlled — from  linear  polar¬ 
ization,  through  elliptical,  to  circularly  polarized  high  har¬ 
monics — without  compromising  the  efficiency  of  the 
upcon version  process  [45,  46].  When  high  harmonics  are 
driven  by  a  bi-chromatic  laser  field  in  which  the  two  spectral 
components  (e.g.  red  (800  nm)  and  blue  (400  nm))  are  cir¬ 
cularly  polarized  with  opposite  helicity  [27,  44,  47],  the 
resulting  HHG  spectrum  consists  of  pairs  of  circularly 
polarized  harmonics — also  with  opposite  helicity.  The  exact 
spectral  position  and  separation  of  the  harmonic  orders 
depend  on  the  driving  laser  wavelengths:  when  the  funda¬ 
mental  and  second  harmonic  of  Ti: sapphire  are  used,  some 
harmonics  rotate  with  the  same  helicity  as  the  red  pump 
(orders  3m  +  1  with  m  =  1,  2,  3  ...)  while  other  harmonics 
rotate  with  the  same  helicity  as  the  blue  pump  (orders 
3m  —  1).  Harmonics  3m  are  forbidden  [27,  29].  The  HHG 
spectrum  in  [46]  is  non-chiral  although  it  is  comprised  of 
circularly  polarized  harmonics.  The  average  chirality  cancels 
since  pairs  of  counter-rotating  circularly  polarized  harmonics 
have  comparable  power.  In  the  time-domain,  the  HHG  fields 
that  correspond  to  such  spectra  consist  of  linearly  polarized 
attosecond  pulse  trains,  with  rotating  polarization  axis 
[30,  46,  48].  Following  the  experimental  generation  of  cir¬ 
cularly  polarized  high  harmonics  [45,  46]  phase  matching  and 
downstream  application  of  circularly  polarized  HHG  was 
demonstrated  in  the  extreme  UV  (EUV)  [47]  and  soft  x-ray 
region  [48],  the  associated  motion  of  the  electrons  during  their 


production  was  explored  theoretically  and  experimentally 
[49],  and  a  novel  3D  attosecond  metrology  technique  for 
characterizing  the  time-domain  field  associated  with  circularly 
polarized  high  harmonics  was  developed  and  demonstrated 
experimentally  [50].  Other  recent  work  generated  circularly 
polarized  high  harmonics  using  crossed  laser  beams,  to  con¬ 
trol  the  propagation  direction  of  left  and  right  circularly 
polarized  HHG  beams  [51].  Finally,  other  techniques  were 
experimentally  demonstrated  to  yield  generation  of  useful 
elliptically  polarized  high  harmonics  [52,  53]. 

An  interesting  question  to  consider  is  the  degree  of  cir¬ 
cular  polarization  of  this  novel  light  source  and  what  factors 
determine  it.  In  [47],  we  measured  nearly  ~100%  degree  of 
circularity  for  the  3m  +  1  harmonics.  Importantly,  we  found 
that  phase  matching  is  helicity  dependent  and  as  such,  it  may 
be  used  for  generating  helically  polarized  attosecond  pulses 
(it  was  recently  discovered  that  the  single  atom  physics  of 
circularly  polarized  HHG  from  p-state  atoms  is  also  helicity 
dependent  [54-56]).  In  that  work,  we  derived  the  phase 
matching  condition  of  the  process  by  assuming  three  con¬ 
servation  laws  for  the  HHG  and  pump  photons:  conservation 
of  energy,  linear  momentum  and  spin  angular  momentum. 
The  same  method  was  used  recently  for  phase  matching  of 
circularly  polarized  x-ray  HHG  [48].  Moreover,  [48]  con¬ 
sidered  the  full  propagation  of  the  HHG  and  driving  laser 
fields,  while  also  including  the  influence  of  the  degree  of 
polarization  of  the  driving  laser  fields,  to  predict  the  HHG 
pulse  duration,  detailed  spectrum  and  hellicity. 

In  this  invited  paper,  we  explore  phase  matching  of  cir¬ 
cularly  polarized  HHG  driven  by  circularly  polarized  counter¬ 
rotating  bi-chromatic  pumps  using  two  analytical  approaches. 
First,  we  extend  our  previous  analytical  approaches  which 
assumed  conservation  of  energy,  linear  momentum  and 
angular  momentum  in  the  exchange  process  between  the 
pump  (fundamental  and  its  second  harmonic)  photons  and  the 
harmonic  photons  [47],  to  general  bi-chromatic,  circularly 
polarized,  counter-rotating  pumps  [48].  Second,  we  derive  a 
more  general  phase  matching  condition  by  requiring  a  new 
propagation  matching  condition  between  the  classical  vec¬ 
torial  field  of  the  bi-chromatic  pump  and  the  resulting  circu¬ 
larly  polarized  high  harmonics.  This  approach  allows  us  to 
identify  that  phase  matching  couples  conservation  of  energy, 
linear  and  spin  momenta.  Under  phase  matching  conditions, 
any  two  (out  of  the  three)  conservation  laws  inevitably  yield 
the  third  one.  We  use  the  classical  fields  approach  to  explore 
the  influence  of  the  pulse  envelopes  on  phase  matching — a 
feature  that  cannot  be  captured  within  the  first  approach. 
Specifically,  we  derive  conditions  for  full-  and  quasi-phase 
matching  of  circularly  polarized  HHG  driven  by  a  long  cir¬ 
cularly  polarized  pulse  and  a  short  counter-rotating  pulse. 
Finally,  we  discuss  the  generation  of  phase-matched  high 
harmonics  that  can  support  helically  polarized  attosecond 
pulses.  We  also  present  experimental  spectra  and  show  that 
they  can  support  attosecond  pulses  with  duration  as  short  as 
120  attoseconds  per  pulse,  and  degree  of  circular  polarization 
as  high  as  80%. 

The  paper  is  organized  as  follows:  in  section  2  we  discuss 
phase  matching  of  high  harmonics  driven  by  bi-chromatic 
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Figure  1.  Experimental  spectrogram  of  HHG  pumped  by  a  bi-chromatic  field  versus  the  relative  phase  between  its  two  spectral  components, 
(a)  For  a  co-linear  bi-chromatic  pump  field  EBC(t)  oc  x  9te[exp(i ut  +  i<^)+exp(2i ujt  +  i</>2u;)]  the  HHG  spectrum  oscillates  with  the 
relative  phase,  <f)2u  —  cj)u,  since  the  shape  of  the  field,  EBC(t)  is  phase  dependent  (see  inset:  red  curve  for  0  phase  and  blue  curve  for  7t/2 
phase),  (b)  HHG  driven  by  counter-rotating  circularly  polarized  bi-chromatic  field,  EBC  ( t )  oc  SKe  [eL  exp  (icjt  +  i <j)u)  +  eR  exp 
(2iojt  +  i</>2u;)]’  is  independent  of  the  relative  phase,  </>2cj  —  <f>u,  where  eL  and  eR  are  left-  and  right-rotating  circularly  polarized  polarizations, 
respectively.  When  varying  the  relative  phase,  the  three-fold  rosette  shape  of  such  a  bi-chromatic  field  (see  inset  for  0  and  n/2  phase) 
undergoes  a  mere  rotation,  which  is  spectrally  undetectable.  We  associate  the  field  shape  with  a  Lissajous  curve  since  the  field  polarization  is 
planar,  rather  than  a  line. 


pump  laser  fields,  focusing  on  the  difference  between  circu¬ 
larly  and  non-circularly  polarized  pump  fields.  Section  3 
offers  a  concise  derivation  of  phase  matching  of  circularly 
polarized  harmonics  using  conservation  laws  for  photons,  as 
derived  first  in  [47].  In  section  4  we  derive  the  phase 
matching  condition  when  treating  the  pump  and  HHG  as 
classical  quasi-monochromatic  vectorial  fields.  We  utilize  the 
derivation  of  section  4  to  investigate  previously  inaccessible 
regimes.  As  an  example,  we  analyze  phase  matching  regimes 
of  circularly  polarized  HHG  pumped  by  a  bi-chromatic  field 
comprised  of  a  short  pulse  and  a  quasi-monochromatic  field: 
the  conditions  for  full  phase  matching  and  quasi-phase 
matching  are  analyzed  in  sections  5  and  6,  respectively. 
Section  7  presents  experimental  highly  chiral  HHG  spectrum 
and  discusses  the  attosecond  pulses  that  this  spectrum  can 
correspond  to. 


2.  Polarization  in  full  phase  matching  of  HHG  by  bi- 
chromatic  drivers 

Phase  matching  of  HHG  driven  by  circularly-polarized  bi- 
chromatic  laser  fields  is  very  different  from  the  case  of  line¬ 
arly  polarized  lasers  [47,  48].  In  HHG  driven  by  bi-chromatic 
linearly  polarized  light,  the  relative  phase  between  the  two 
spectral  components  of  the  driving  field  dramatically  influ¬ 
ences  the  phase  and  the  intensity  of  the  HHG  field  [57-59], 
because  it  changes  the  maximal  instantaneous  field  and  the 
shape  of  the  bi-chromatic  driving  field.  As  a  result,  the  rela¬ 
tive  phase  modifies  the  tunnel  ionization  rate  of  the  atom  [60], 


the  trajectories  of  these  electrons  in  the  continuum,  and  the 
resulting  HHG  yield  [61-63].  Figure  1(a)  presents  the  HHG 
spectrum  as  a  function  of  the  relative  phase  between  the  two 
colors  of  a  co-linearly  polarized  bi-chromatic  driving  laser 
focused  into  a  gas  jet  (see  apparatus  in  [46]).  The  inset 
illustrates  the  temporal  shape  of  the  field  for  relative  phase  of 
0  and  7t/2  between  the  two  spectral  components.  When  a  co- 
linearly  polarized  bi-chromatic  field  propagates  in  a  dis¬ 
persive  medium  the  phase  of  the  two  spectral  components 
travel  at  different  speeds,  so  their  relative  phase  cannot  be 
fixed.  The  amplitude  and  phase  of  the  resulting  HHG  field 
oscillate  along  the  propagation  direction  with  a  periodicity 
determined  by  the  refractive  indices  of  the  two  spectral 
components.  The  consequence  of  these  oscillations  is  that  full 
phase  matching  of  high  order  harmonics  is  extremely  chal¬ 
lenging  [64]. 

On  the  other  hand,  in  the  case  of  a  bi-chromatic  circularly 
polarized  driving  laser  field,  the  phase  between  the  two 
spectral  components  merely  rotates  the  rosette-shaped  field  in 
the  polarization  plane  (see  inset  of  figure  1(b)).  The  relative 
phase  does  not  change  the  peak  intensity,  nor  the  temporal 
shape  of  a  pump  field  and  not  the  resulting  HHG  spectrum. 
Figure  1(b)  presents  a  spectrogram  of  HHG  driven  by  a 
counter-rotating  bi-chromatic  pump  field  (see  apparatus  in 
[46]).  Clearly,  the  HHG  spectrum  is  independent  of  the 
relative  phase  between  the  pump  spectral  components.  In  this 
respect,  phase  matching  of  circularly  polarized  HHG  driven 
by  a  circularly  polarized  pump  is  similar  to  phase  matching  of 
circularly  polarized  harmonic  generation  driven  by  a  linearly 
polarized  laser  rotating  in  an  optically  active  medium 


3 


J.  Phys.  B:  At.  Mol.  Opt.  Phys.  49  (2016)  123501 


Invited  Article 


b) 


kq  =  +  1  j  k1  +  mk2 


o 


— ► — ►-y 

k±  kt  k±  k2 


k. 


kn  =  (m  —  1\  ki  +  mk2 

O  ' 


-► - ► — 0- 


Photon  linear  momentum  (up  to  h ) 


-► 


A  k 


Figure  2.  Conservation  of  energy,  spin  angular  momentum  and  linear  momentum  in  the  upcon version  of  many  circularly  polarized  u\  and  uj2 
photons  into  a  HHG  photon  in  the  qth  harmonic  order,  ujq.  (a)  HHG  orders  uq  =  (m  ±  \)u\  +  mu 2  conserve  odd  parity,  spin-angular 
momentum  and  energy.  The  emitted  HHG  photon  inherits  the  spin  of  the  extra  pump  photon  (extra  u\  marked  red  and  extra  u2  marked  blue), 
(b)  The  upconversion  process  in  a  dispersive  medium  is  helicity  selective — if  one  helicity  is  phase  matched,  the  phase  of  the  opposite  helicity 
is  mismatched  (marked  by  A k). 


[33,  65].  In  both  cases,  the  driving  field  rotates  during  pro¬ 
pagation,  yet  its  shape  and  amplitude  are  conserved.  Notably, 
the  fact  that  the  phase  of  circular  HHG  is  relatively  insensitive 
to  intensity  variations  (compared  to  linear  HHG)  [30]  makes 
it  more  robust  to  pulse  envelop  effects  and  beam  focusing 
dynamics. 

3.  Full  phase  matching  condition  from  photon 
conservation  laws 

Below,  we  generalize  the  derivation  of  the  phase  matching 
equation  in  [47]  to  the  case  of  any  bi-chromatic  pump;  not 
necessarily  a  fundamental  and  its  second  harmonic,  as  was 
recently  experimentally  demonstrated  for  soft  x-ray  HHG 
driven  by  1300  nm  and  800  nm  circularly -polarized  lasers 
[48].  This  section  addresses  the  case  where  the  pulses  are 
long,  so  temporal  variations  in  their  envelope  can  be 
neglected.  When  HHG  is  driven  by  a  bi-chromatic  field 
composed  of  co-propagating  left-  and  right-rotating  circularly 
polarized  fields,  of  frequency  uj\  and  cu2,  respectively,  the 
photons  conservation  laws  are  given  by: 

energy  uq  =  £iu\  +  mu>2,  (1) 

spin  angular  momentum  i  =  m  db  1 ,  (2) 

linear  momentum  kq  —  ik\  +  mk2,  (3) 

t  and  m  are  the  number  of  uj\  and  u2  photons  annihilated, 
respectively,  from  the  bi-chromatic  pump  to  generate  a  HHG 
photon  at  the  harmonic  order  q  with  frequency  uoq  —  qu\.  If 
u2/u\  is  an  integer,  so  is  q.  The  HHG  photon,  uq,  inherits  the 
spin  of  the  one  extra  photon  it  acquired,  and  is  therefore,  left- 


rotating  for  £  =  m  +  1  and  right-rotating  for  £  =  m  —  1 .  The 
conservation  of  energy  and  spin  angular  momentum  are 
illustrated  in  figure  2(a).  The  momentum  of  an  HHG  photon, 
kq,  (the  factor  of  H  is  ommited)  originates  from  annihilation  of 
l  photons  with  momentum  k\  and  m  photons  with  momentum 
k2  (equation  (3)).  Conservation  of  linear  momentum  is  the 
condition  for  phase  matching  of  circularly  polarized  harmo¬ 
nics;  the  HHG  process  is  fully  phase  matched  when 
equation  (3)  is  fulfilled.  However,  due  to  chromatic  and 
plasma  dispersion  in  the  generating  medium,  the  fields  at  uj i 
and  u2  have  different  refractive  indices  so  k2  /  k\^  uj2/uj\. 
The  momentum  mismatch  (also  called  phase  mismatch),  A k, 
is  retrieved  from  equations  (1),  (2)  and  (3)  to  be 

A  k  =  £k\  +  mk2  —  kq  =  £—Ani  +  m—An2.  (4) 

c  c 

Here  A k  is  the  phase  mismatch  of  HHG  with  frequency 
uq  =  £uj\  +  mu2,  where  t  =  m  =L  1.  The  refraction  terms 
Ani 2  are  given  by  A n^2  =  (^i;2  —  nq),  where  nlj2  is  the  index 
of  refraction  including  plasma,  modal  or  Gouy  phase  terms  of 
the  pumps  and  uj2,  respectively  [66,  67],  and  nq  is  the  index 
of  refraction  for  the  harmonic  order  q.  The  condition  for  full 
phase  matching  is: 

cj2A n2  _  l  _  (mil) 
u\Ani  m  m 

According  to  equation  (5),  A n2  and  An\  should  have 
opposite  signs.  Equations  (3)-(5)  lead  to  several  conclusions: 
First,  only  one  HHG  order  can  be  fully  phase  matched  for 
given  dispersion  terms  Ani^2.  Let  us  assume  that  the  phase 
matched  harmonic  rotates  like  the  field.  Other  co-rotating 
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HHG  orders  (i.e.  orders  t  =  m  +  1)  differ  by  a  small 
momentum  addition  c~l(uJiAni  +  cj2An2),  where  the  phase 
mismatch  grows  for  HHG  orders  further  away  from  the 
fully  phase  matched  HHG  order.  Note  that  the  term 
(<jj\  Any  +  ^2  An2)  is  small  since  An\  and  A n2  have  opposite 
signs.  Third,  all  counter-rotating  HHG  orders  (£  =  m  —  1) 
have  a  larger  phase  mismatch  than  the  co-rotating  HHG 
orders.  When  compared  with  the  phase  matched  HHG  order, 
HHG  orders  with  opposite  helicity  lack  two  uj\  photons  so 
their  phase  mismatch  is  approximately  Ak  «  2c~lu\An\. 
Indeed,  usually  uj\An\  (u\An\  +  uo2An2)  so  co-rotating 
harmonics  are  better  phase  matched  than  counter-rotating 
orders.  Figure  2(b)  illustrates  the  helicity  selectiveness  of 
such  a  phase  matching  scheme  where  left-rotating  HHG  is 
phase  matched  and  right-rotating  HHG  exhibits  a  large  phase 
mismatch  (see  also  figure  1  in  [47]). 

The  use  of  conservation  laws  in  the  above  derivation 
implicitly  utilizes  continuous  symmetries  of  the  system  in  an 
intuitive  way.  The  reliance  of  the  phase  matching  equation  on 
conservation  of  energy,  linear  momentum  and  angular 
momentum  indicates  that  the  system  (light  and  isotropic 
matter)  is  invariant  to  temporal  delays,  spatial  translations  and 
rotations,  respectively.  However,  one  can  address  the  sys¬ 
tem’s  symmetries  directly.  In  the  next  section,  we  derive  the 
phase  matching  equation  using  symmetries  of  the  classical 
pump  field  propagating  in  an  isotropic  medium  and  retrieve 
the  same  phase  matching  condition  as  in  this  section,  i.e. 
equation  (4). 


4.  Full  phase  matching  condition  by  vectorial  fields 
symmetries. 

In  this  section,  we  derive  the  condition  for  full  phase 
matching  by  analyzing  the  operator  that  propagates  the  vec¬ 
torial  pump  field  in  circularly  polarized  HHG  (we  assume 
nothing  about  the  details  of  HHG  or  about  the  photons  that 
participate  in  the  process).  As  will  be  shown  in  this  and  in 
following  sections,  this  formalism  is  much  more  general  than 
the  photon  conservation  laws  approach.  To  derive  the  con¬ 
dition  for  phase  matching,  we  identify  a  unitary  propagator 
for  the  bi-chromatic  pump  in  a  chromatically  dispersive 
medium  and  find  the  condition  for  which  the  emitted  HHG 
field  builds  up  coherently  throughout  the  entire  generation 
medium.  Notably,  unlike  phase  matching  of  linearly  polarized 
HHG,  where  pump  and  HHG  fields  have  the  same  polariza¬ 
tion,  here  the  fields  cannot  be  represented  by  complex-valued 
scalars.  Our  description  accounts  for  the  vectorial  nature  of 
the  fields. 

We  start  by  assuming  a  bi-chromatic  field,  composed  of 
circularly  polarized,  counter-rotating  components  at  fre¬ 
quencies  uj\  and  cj2  that  propagates  along  z  direction  in  an 
isotropic  medium.  For  simplicity,  we  neglect  diffraction  and 
solve  for  time  and  propagation  axes  and  assume  plane  wave 
propagation.  We  define  z  =  0  in  an  arbitrary  position.  The 


bi-chromatic  field,  EBC,  at  z  =  0  is  given  by: 

gBc(T,z  =  0)  =  ^(j)e-  +  ^(_'i)el-. 

(6) 


Here  E\  and  E2  are  scalar  complex  amplitudes  of  the  two 
spectral  components  of  the  bi-chromatic  field.  The  2x1 
vector  is  the  projection  of  the  electric  field  on  axes  v  and  y  in 
the  polarization  plane,  r  represents  time  shifts  relative  to  a 
mathematical  frame  of  reference  traveling  at  the  phase 
velocity  of  cjq,  vq  =  c/nq ,  where  c  is  the  speed  of  light,  so 
r  =  t  —  nqz/c.  This  becomes  useful  when  we  look  for  a 
coherent  buildup  of  harmonic  order  q — by  using  r  instead  of 
t,  the  phase  of  a  traveling  HHG  field  is  z-independent.  For  a 
propagating  bi-chromatic  field,  the  phase  term, 
takes  the  form  Q^i,2r~ikh2z- Anh2  when  we  use  r,  where 
k\  2  =  u\ ,2/c  and  A n1>2  =  n\,2  —  nq.  Therefore,  the  bi-chro- 
matic  field  at  z  >  0  is 


EBC  (t,  z)  =  ffte 

|  1  \gioo2T— i^-zAn2 

V2  V— i/ 


(7) 


Note  that  the  dispersion  terms  n\  2  may  include  nonlinear 
contributions  (e.g.  Kerr  and  self-generated  plasma)  as  long  as  the 
field  intensity  is  constant  throughout  propagation.  It  is  instructive 
to  obtain  the  propagation  operator  that  transfer  £BC  (r,  0)  of 
equation  (6)  to  EBC  (r,  z)  of  equation  (7).  To  do  so,  we  rearrange 
the  exponential  arguments  for  uj\  and  uj2,  as  uj\  (r  +  t)  +  6 
and  uj2  (t  +  t)  —  <5,  respectively,  and  find  the  proper  t,  6 


Hz) 

Hz) 


z  (viAni  +  cj2Au2) 


C  {id\  ~\~  U 2) 

- (ni  - 

c  (jj\  T  uj2 


ni). 


(8) 

(9) 


Using  t  and  6  (equations  (8)  and  (9))  in  equation  (7)  leads 
to: 

Esc ( r ,  z)  =  eHz)frR6(z)EBC (r,  0).  (10) 

Here  e^  is  a  unitary  time  delay  operator  and  R$  is  the  unitary 
rotation  matrix  in  the  polarization  plane  given  by 


(  cos  6  (z)  sin  6  (z)  \ 
\-sin  8(z)  cos  8(z)\ 


(ID 


The  time  delay,  t ,  can  be  thought  of  as  resulting  from  an 
effective  refractive  index  of  the  bi-chromatic  field, 
t  =  —  (z/c)A«effeCtive,  where  Arc|f(ictive  is  a  weighted  average 
of  An\  2  with  u i>2  as  weights. 

In  HHG,  a  process  is  fully  phase  matched  if  emissions 
from  any  two  propagation  distances  interfere  constructively, 
regardless  of  the  propagation  distance,  Az,  between  them. 
Since  the  bi-chromatic  pump  field  driving  the  HHG  emission 
merely  rotates  and  delays  as  it  propagates,  the  HHG  fields 
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Figure  3.  Propagation  of  a  ‘rosette  Lissajous  curve’  of  a  bi-chromatic  pump  field  with  frequencies  u\  and  uj2  =  2cjh  Due  to  a  difference  in  the 
indices  of  refraction  for  uo x  and  uj2,  the  rosette  shape  rotates  as  it  propagates.  Insets  (b)  and  (c)  show  the  difference  in  the  rosette  rotation  and 
time  delay  at  two  position  on  the  propagation  axis,  separated  by  A z.  The  black  arrows  in  the  insets  is  the  electric  field  at  r  =  0.  t  corresponds 
to  the  delay  accumulated  over  Az  between  the  rosette  and  HHG  fields  due  to  their  different  phase  velocities.  Circularly  polarized  HHG  are 
fully  phase  matched  when  the  phase  differences  induced  by  the  time  delay  and  by  the  rotation  of  the  bi-chromatic  pump,  t  and  ^rotation, 
respectively,  exactly  counteract  each  other.  Under  such  a  balance,  the  HHG  emitted  throughout  the  propagation  builds  up  coherently.  The 
screw-like  propagation  of  the  bi-chromatic  field  is  a  visual  demonstration  for  the  break  of  chiral  symmetry,  which  is  required  for  helicity 
selective  HHG. 


emitted  at  different  positions  are  identical  up  to  a  rotation 
angle  (same  rotation  angle  as  the  driver  field)  and  time  delay. 
Thus,  for  full  phase  matching  of  g-order  harmonic  field 
generation,  the  net  effect  of  the  rotation  and  time  delay  should 
yield  the  identity  operation  (see  figure  3): 

Eqc^T  =  (12) 

Conveniently,  the  use  of  r  instead  of  t  eliminates  any 
explicit  z-dependency  of  the  HHG  field.  Equation  (12)  is  a 
simple  eigenvalue  equation.  Its  solutions  are  the  polarization 

eigenvectors  Eq  =  that  correspond  to  left-rotating  (+) 

and  right-rotating  (— )  circular  polarization,  with  eigenvalues 
eiu^?±i<§,  respectively.  Therefore,  the  phase  matching  equation 
condition  for  ujq  is 

uqt  ±8  =  0,  (13) 

where  t  and  8  are  defined  in  equations  (8)  and  (9).  The  phase 
matching  equation  (13)  can  be  written  explicitly  as 

0  =  - - - [0Oq  (Q0\Ari\  ±  CJ2An2) 

c(  LJ\  +  CJ2) 

±  ujiu2{/±n\  —  An2)].  (14) 

Equations  (12)— (14)  correspond  to  the  case  of  full  phase 
matching.  It  is  also  instructive  to  formulate  the  phase- 


mismatch  situation.  In  this  case,  we  add  a  phase 
mismatch  term,  e_lMAz,  to  equation  (12)  and  obtain: 
EqeluJ<iT-lAkAz  =  ±  ^Az)  drR^^EqO}^ .  Hence,  the  phase 
mismatch  generalization  to  equation  (14)  is  given  by: 

7 r  1 

—  =  A k  = - [uq(uJiAni  ±  u2An2) 

Lc  c(lu\  +  cj2) 

db  LJiCJ2(Ani  —  A n2)],  (15) 

where  Lc  is  the  coherence  length  of  the  conversion  process. 

Thus  far  we  derived  two  different  equations  for  phase 
matching  of  circularly  polarized  HHG  driven  by  bi-chromatic 
circularly  polarized  counter-rotating  drivers:  equations  (4) 
and  (15).  Equation  (15)  is  more  general  as  it  does  not  require 
to  ‘count’  the  number  of  annihilated  photons,  l  and  m,  as 
input  parameters.  The  two  equations  coalesce  if  one  assumes 
any  two  out  of  the  three  conservation  laws  (equations  (1),  (2) 
or  (3)).  For  example,  equation  (3)  can  be  derived  directly 
from  equations  (15),  (1)  and  (2).  In  a  similar  fashion, 
equation  (2)  can  be  derived  from  equations  (15),  (1)  and  (3). 
This  feature  reveals  that  the  phase  matching  in  our  scheme 
couples  the  conservation  laws  for  the  photons  in  the  process. 

Figure  3  illustrates  the  propagation  effects  on  a  circularly 
polarized  counter-rotating  bi-chromatic  field  comprising  a  fun¬ 
damental  field,  u i,  and  its  second  harmonic,  uo2  =  2uj\.  As  the 
field  travels  in  a  chromatically  dispersive  medium,  the 
three-fold  rosette  shape  of  the  field  rotates  since  uj\  and  uj2 
acquire  a  phase  shift  (figure  3(a)).  This  screw-like  propagation 
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clearly  shows  that  the  pump  is  chiral,  and  therefore,  one  should 
expect  a  major  difference  in  the  coherence  lengths  of  high 
harmonics  with  left  helicity  versus  right  helicity.  Circularly 
polarized  fields  are  a  unique  outcome  for  such  a  screw-like 
propagation  since  spatial  rotation  simply  induce  a  phase 
[33,  49,  65].  According  to  equation  (13),  the  HHG  phase  is 
matched  when  the  chirally  selective  phase  of  the  spatial  rota¬ 
tion,  <5,  balances  the  time  delay  phase,  the  result  should  be  ujqt, 
which  is  a  non-chiral  effect  of  the  propagation.  Figures  3(b)  and 
(c)  show  6  and  t  for  two  positions  along  the  propagation  axis, 
z  =  0  and  z  =  Az.  Notably,  while  figure  3  illustrates  the  bal¬ 
ance  between  the  time  delay  and  rotation  for  the  case  of 
uo2  =  2aq,  it  works  for  any  ^  cj2. 


5.  Full  phase  matching  of  circularly  polarized  HHG 
by  a  short  and  a  long  pump  pulses 


In  the  previous  sections,  we  derived  the  phase  matching 
equation  while  neglecting  the  envelopes  of  the  bi-chromatic 
pump  pulses.  In  this  section  we  extend  the  analysis  to  include 
the  effect  of  pulse  envelope.  Specifically,  we  shall  discuss  full 
phase  matching  conditions  for  circularly  polarized  HHG 
pumped  by  a  combination  of  a  short  circularly  polarized  pulse 
at  uj\  and  a  long  counter-rotating  pulse  at  uj2  [54].  The 
properties  of  the  HHG  formed  by  such  a  pump  field  depend 
strongly  on  the  carrier-envelope  phase,  that  is,  the  time  dif¬ 
ference  between  the  pulse  envelope  and  the  peak  of  the 
oscillating  field.  Models  based  on  continuous  symmetries  or 
conservation  laws  for  photons  (e.g.  the  model  in  section  3) 
neglect  this  dependence  and  are  therefore  incomplete.  The 
vectorial  fields  model  (section  4),  on  the  other  hand,  can  be 
extended  to  include  the  effect  of  the  pulse  envelope. 

In  short,  if  the  velocities  of  the  envelope  and  the  peak  of 
the  carrier  bi-chromatic  field  are  equal  (although  the  carrier  field 
might  rotate),  then  the  shape  of  the  field  is  independent  of  the 
propagation.  In  this  case,  high  harmonics  that  propagate  with 
the  same  velocity  would  be  fully  phase  matched.  If,  on  the  other 
hand,  the  velocities  of  the  envelope  and  carrier  field  differ,  then 
the  field’s  shape  oscillates  as  it  propagates  through  the  medium, 
and  the  HHG  field  may  only  be  quasi  phase  matched  (see 
section  6).  In  this  and  the  following  section  the  vectorial  nature 
of  the  pump  field  is  addressed  explicitly  and  accounts  for  var¬ 
iations  of  orientation,  amplitude  and  phase  of  the  electric  field 
throughout  propagation.  For  simplicity,  we  take  into  account 
the  group  and  phase  velocities  of  the  pump  fields  but  leave  out 
higher  order  effects  such  as  group  velocity  dispersion,  and 
effects  of  ionization  within  the  duration  of  the  short  pump  pulse. 

We  assume  that  the  bi-chromatic  field,  E^c,  entering  the 
medium  at  z  =  0  is  given  by: 


e£c  (r,  0)  =  Vie 


(16) 


The  field  at  frequency  uj\  is  confined  to  a  Gaussian 
envelope  with  a  standard  deviation  at,  and  the  pulse  duration 


of  the  field  at  uo2  is  much  longer  and  is  taken  as  infinite.  Here 
we  marked  the  bi-chromatic  field  with  a  superscript  ‘p’  to 
denote  the  pulsed  nature  of  the  field  and  to  distinguish  it  from 
the  bi-chromatic  fields  in  previous  sections.  The  pulse 
envelop  is  taken  as  Gaussian  for  simplicity,  however,  the 
derivation  is  general  and  can  be  applied  to  any  pulse  envel¬ 
ope.  Due  to  propagation  of  distance  z  in  the  medium,  the  field 
undergoes  chromatic  dispersion  (similar  to  equation  (7)) 
while  the  pulse  envelope  aquires  a  group  delay  z/vg,  so 


Ebc  (r>  z)  =  91 e 


-^-Qe“,rw“iZ  An'~  '  i)) 


,iu)2T—  lk2Z  A«2 


vg  is  the  group  velocity  of  a  pulse  centered  at  u)\  and  we 
substituted  r  =  t  —  z/vq,  and  An\^2  =  (n\2  —  nq).  With  the 
same  logic  as  in  section  4,  we  wish  to  extract  the  symmetry 
propagator  that  links  E^c  (r,  0)  and  E^c  (r,  z).  We  claim  that 
the  HHG  phase  is  matched  when  the  field  of  the  HHG  order  q 
conforms  to  the  same  propagator.  We  start  by  extracting  the 
time  and  rotation  operators,  as  in  equation  (10) 


Ebc  (r>  z) 


By  using  t  and  6  we  eliminate  the  terms  e_1^2Z  Ani’2  but 
add  an  explicit  time  shift,  7,  to  the  pulse  envelope.  For  the 
operator  to  be  unitary,  as  in  equation  (12),  equation  (13) 

should  hold  and  also  the  term  t  +  z\  — - -I  should  nullify, 

\vg  vqJ 


leading  to  the  condition 


0  =  t  +  z 


1 


1 


V(U 


yq 


=  z 


1  (cjiA^i  +  u2An2) 

C  ((jJi  +  lj2) 

(19) 


This  additional  condition  for  full  phase  matching  cannot 
be  predicted  by  considering  the  conservation  of  energy,  spin 
angular  momentum  and  linear  momentum  of  the  pump  and 
HHG  photons,  and  requires  the  explicit  approach  presented 
here.  Equation  (19)  accounts  for  the  relative  time  delay 
between  the  bi-chromatic  field  (induced  by  chromatic  dis¬ 
persion)  and  the  group  delay  of  the  pulse  envelope.  When 
equation  (19)  is  fulfilled  and  this  relative  time  delay  nullifies, 
the  shape  of  the  bi-chromatic  field  is  independent  throughout 
propagation,  up  to  trivial  time  delays  and  rotations.  The  HHG 
field  that  this  field  produces  will  be  identical  along  propa¬ 
gation,  up  to  these  delays  and  rotation,  similarly  to  the  case 
described  in  section  4.  Therefore,  the  phase  matching  con¬ 
dition  for  the  HHG  field  is  equal  to  the  condition  in 
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equation  (15),  with  the  addition  of  equation  (19) 


Ak 


0 


1 


-K  (wiAni  +  002 An2) 


c(oJi  +  102) 

=b  CjOiCJ2  (A/2i  —  Aft2)] 


(ooiA^i  +  oo2A  n2) 

(lOi  +  102) 


+  c 


1 


(20) 


Although  equation  (20)  is  similar  to  equation  (15),  in  a 
normally  dispersive  medium  equation  (20)  allows  phase 
matching  only  for  harmonic  orders  with  the  same  helicity  as 
uj\.  In  a  normally  dispersive  medium  n2  >  n\,  vq  «  c  and 
c  >  vg,  so  the  =b  sign  must  be  +  in  order  to  nullify  Ak. 
Intuitively,  full  phase  matching  occurs  only  if  the  envelope 
aquires  the  same  time  delay  as  the  rosette  of  the  bi-chromatic 
field,  and  therefore  the  entire  field  rotates  and  delayes  in  time, 
without  changing  its  shape.  Note  that  these  are  indeed  the 
dominant  HHG  orders  in  the  chiral  spectra  that  we  measured 
in  [47]  (figure  (4)  and  also  in  figure  (4)  of  this  work).  In  the 
next  section  we  show  that  if  the  shape  of  the  pumping  bi- 
chromatic  field  does  vary  with  propagation,  then  the  HHG 
process  can  be  quasi  phase  matched  (QPM). 


6.  Quasi  phase  matching  of  circularly  polarized  HHG 

In  this  section  we  consider  the  case  where  the  envelope  and 
the  rosette  field  do  not  propagate  with  equal  velocities  (i.e. 
equation  (19)  is  not  satisfied),  leading  to  a  continuous  slip 
between  the  phase  and  the  envelope  of  the  bi-chromatic 
pump.  In  such  a  case,  the  phase  and  amplitude  of  the  emitted 
HHG  field  vary  periodically  as  the  pump  propagates  in  the 
generation  medium  and  HHG  can  be  QPM  [67],  but  not  fully 
phase  matched.  To  find  the  periodicity  of  the  phase  slip  we 
use  a  property  of  counter-rotating  bi-chromatic  pump:  a  time 
delay  of  2tt/(cji  +  uj2)  simply  rotates  the  field  by  a  spatial 
angle  2ttu;i/(uji  +  uo2).  A  rotation  by  this  angle  does  not 
change  the  bi-chromatic  field.  In  the  particular  case  of 
cj2  =  2uq,  a  delay  of  27r/3o;i  (1/3  of  an  optical  cycle)  is 
equivalent  to  a  rotation  by  120°,  which  is  exactly  the  rotation 
symmetry  of  the  bi-chromatic  field.  Therefore,  although  the 
envelope  delay  term  in  equation  (18),  t  +  z(vg_1  -vql),  may 
not  nullify  continuously,  it  nullifies  after  propagating  one 
period,  marked  A  (see  notation  in  section  2.9  in  Boyd  [68]). 
For  example,  if  the  carrier  bi-chromatic  field  and  the  envelope 
reach  their  maxima  together  at  some  position,  they  also  reach 
their  maxima  jointly  after  propagating  A,  albeit  with  some 
time  delay  and  rotation.  It  is  convenient  to  define  a  time 
delay,  ta,  and  to  write  its  relation  with  A  in  the  following 
form 


derivation,  since  it  takes  into  accout  the  temporal  and 
rotational  symmetries  of  the  the  bi-chromatic  field. 

We  write  the  bi-chromatic  field  after  it  propagates  A 
using  equation  (18).  We  employ  a  unitary  operation  that  adds 
ta  within  the  brackets  and  adds  —  ta  to  the  time  delay  operator 
to  simplify  the  envelope  group  delay  term  (see  equation  (21)) 

E\  / 1\  iwiCr+TA)--^ 

vr(i)e 

+  //(  jeiU2(r+TA)l  (22) 

The  additional  delay,  ta,  is  apparent  in  the  phasors 
eiu;u(T+TA)  xhe  definition  t\  =  27t/(cji  +  uj2)  now  becomes 
useful,  since  according  to  the  symmetries  of  the  bi-chromatic 
field,  a  spatial  rotation  by  2ttuji/(uji  +  cj2)  removes  ta  from 
the  phasors  and  retreives  E^c(r,  0)  in  the  square  brackets. 
Thus,  we  can  write  a  unitary  discrete  propagator  for  the  bi- 
chromatic  field,  including  the  effects  of  the  envelope 

£bpc  (r,  A)  =  (e(?“rA) d^R6+bHT j  £BPC  (r,  0).  (23) 

The  quasi-phase  matching  condition  for  the  HHG  field, 
Eq,  with  frequency,  ujq,  can  be  retreived  with  the  same  logic 
leading  to  equations  (12)  and  (15).  To  be  QPM,  the  HHG 
field  should  be  equal  at  two  positions  separated  by  distance  A, 
where  deviations  from  QPM  are  expressed  by  the  phase 
mismatch  term,  Akq 

EqeitoqT-iAkQA  =  ^-T^R6+UiTKEqe^T.  (24) 

The  condition  for  QPM  is  retreived  by  substituting  cir¬ 
cularly  polarized  HHG  into  equation  (24) 

—  AkgA  +  27 rl  =  (cjqt(A)  zb  6(A))  +  (— =b  gji)ta, 

(25) 

where  l  is  an  integer  representing  the  QPM  order  (not  to  be 
confused  with  l  of  equation  (1)).  We  use  equations  (8),  (9) 
and  (15)  to  replace  the  left  brackets  with  the  phase  mismatch, 
—AkA,  and  use  equation  (21)  to  get  the  explicit  QPM 
equation 

—  AkqA  +  27 r/  =  —AkA  +  27r(—ujq  zb  +  cj2). 

(26) 


£bc  (r,  A)  =  ef-'-^WRsfte 


By  the  relation  uoq  =  +  cj2)  zb  u)\  (equations  (1) 

and  (2))  we  get  a  simple  condition  for  QPM 

2tt 

Ako  =  Ak  +  —  (Z  +  m),  (27) 

A 


27 r 

TA.  =  - 

+  CJ2 


t  (A)  +  A 


(21) 


where  t  is  defined  in  equation  (8).  The  right-hand  side  of  (21) 
allows  us  to  nullify  the  envelope  delay  term  in  equation  (18), 
where  the  definition  ta  =  +  cj2)  is  usful  for  our 


where  /  is  order  of  QPM  and  m  is  related  to  the  HHG  order. 
Equation  (27)  implies  that  for  QPM,  Akg  =  0,  the  phase 
mismatch  from  the  bi-chromatic  carrier  wave,  Ak,  should 
balance  an  integer  order  of  the  periodicity  of  the  carrier 
envelope  phase  slip,  2tt/A.  If  an  harmonic  order  with  a  given 
m  is  quasi-phase  matched  according  to  equation  (27)  with  a 
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Figure  4.  Highly  chiral  HHG  and  attosecond  pulses,  (a)  Chiral  spectrum  recorded  in  gaseous  neon  [47],  where  the  left-rotating  HHG  orders 
(red)  dominate  over  the  right-rotating  HHG  orders  (blue).  Furthermore,  the  HHG  radiation  acquires  additional  chirality  since  (b)  the  degree  of 
circular  polarization,  defined  as  the  normalized  Stokes  parameter  S3  /So,  is  100%  for  the  left-rotating  HHG  and  40%  for  the  right-rotating 
HHG  (see  figure  S.3  in  [47]).  (c)  In  time  domain,  the  above  spectrum  corresponds  to  a  train  of  highly  chiral  attosecond  bursts,  supporting 
pulse  duration  as  short  as  120  attoseconds  (red  line).  Interestingly,  neighboring  attosecond  bursts  have  lower  degree  of  circular  polarization 
(noted  in  %).  Therefore,  if  such  emission  probes  a  chiral  phenomenon  (e.g.  circular  dichroism),  most  of  the  chiral  signal  (blue  line)  comes 
from  a  single  attosecond  pulse,  (d)  Experimentally  acquired  spectrum  from  a  waveguide  full  of  neon  gas  showing  extremely  high  selectivity 
for  left-rotating  harmonics  (red)  over  right-rotating  harmonics  (blue).  The  extremely  high  chirality  of  such  a  spectrum  can  be  directly 
identified  from  the  contrast  between  the  left-  and  right-rotating  HHG  orders.  It  supports  a  train  of  attosecond  pulses  with  similar  features  as 
in  (c). 


given  QPM  order  /,  other  HHG  orders  (different  m,  with  equal 
or  opposite  helicity)  may  still  nullify  AJcq  for  a  different  /, 
however,  such  nullification  is  not  guaranteed. 

To  conclude  this  section,  we  showed  that  quasi  phase 
matching  is  possible,  and  it  is  helicity  selective  (through 
the  ±  term  in  A k)  similarly  to  full  phase  matching.  However, 
it  is  governed  by  a  different  equation.  The  physical  inter¬ 
pretation  of  equations  (25)-(27)  is  more  complex  than  for  the 
case  of  full  phase  matching.  To  understand  this  system 
intuitively,  consider  the  case  where  the  envelope  of  the  field 
at  uj\  overlaps  with  the  peak  of  the  bi-chromatic  field,  and  the 
pulse  is  short  enough  to  include  only  a  single  burst  of  HHG 
radiation.  Since  the  group  velocity  of  the  short  pulse  differs 
from  the  phase  velocity  of  the  field,  after  some  propagation 


distance,  A,  the  pulse  envelope  ovelaps  with  an  adjacent  peak 
of  the  bi-chromatic  field’s  rosette.  At  that  point,  the  shape  of 
the  field  is  as  it  was  originally  up  to  a  temporal  delay,  t  —  ta, 
and  rotation,  uo\T^.  The  HHG  emitted  from  atoms  separated 
by  A  is  identical,  up  to  the  these  delay  and  rotation.  HHG  is 
QPM  if  the  emissions  at  every  A  differ  by  2i rZ  and  interfere 
constructively. 


7.  Helically  polarized  attosecond  pulses 

Generation  of  circularly  polarized  high  harmonics  by  bi- 
chromatic  counter-rotating  circularly  polarized  drivers  results 
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in  a  comb  of  harmonic  pairs,  where  in  each  pair  the  harmonics 
rotate  in  opposite  directions,  and  can  span  the  EUV  and  soft 
x-ray  regions  [46-48].  Generally,  the  average  chirality  of  the 
field  is  rather  small  because  the  intensities  of  the  left-rotating 
and  right-rotating  harmonics  are  comparable  [46].  Thus,  in 
order  to  generate  helically  polarized  attosecond  pulses,  we 
need  to,  somehow,  break  this  symmetry.  We  need  a  helicity- 
dependent  filter.  As  discussed  in  previous  sections,  phase 
matching  of  circularly  polarized  high  harmonics  is  strongly 
helically  dependent  [47]. 

In  this  section,  we  present  experimentally  measured 
spectra  of  circularly  polarized  HHG  with  a  high  degree  of 
average  chirality  and  show  that  this  measured  spectrum  can 
support  attosecond  pulses  with  a  large  degree  of  circular 
polarization  (the  experimental  corroboration  for  phase 
matching  was  presented  in  [47,  48]).  Figure  4(a)  shows  a 
spectrum  for  circularly  polarized  HHG  from  neon  that  appears 
on  figure  3(b)  in  [47].  The  HHG  process  was  driven  by  a  cj\ 
(790  nm)  circularly  polarized  left  rotating  driver  and  its  sec¬ 
ond  harmonic,  i.e.  cj2  =  which  rotated  in  the  opposite 
direction.  In  this  case,  the  allowed  harmonic  orders  are 
q  =  3m  ±  1,  where  m  is  an  integer  and  q  =  uq/uj\.  Har¬ 
monic  orders  q  =  3m  +  1  are  left-rotating,  as  and  HHG 
orders  q  =  3m  —  1  are  right-rotating,  as  uj2,  where  HHG 
orders  q  =  3m  are  suppressed  (marked  red,  blue  and  green, 
respectively,  in  figure  4(a)).  Clearly,  the  left-rotating 
harmonic  peaks  are  3-5  times  higher  than  their  adjacent  right- 
rotating  peaks.  Notably,  we  also  measured  the  degree  of 
circular  polarization  (i.e.  the  normalized  Stokes  parameter, 
S3 /So,  which  is  the  figure  of  merit  for  magnetic  chiral 
experiments.  See  p  345  in  [69])  of  four  harmonics  in  this 
spectrum  (figure  4(b),  taken  from  figure  S.3  in  the  supple¬ 
mentary  information  of  [47]).  As  shown  in  figure  4(b),  the 
degree  of  circular  polarization  of  right-rotating  HHG  orders 
are  ~40%,  whereas  left-rotating  HHG  orders  are  ~100% 
circularly  polarized.  In  short,  the  left-rotating  part  of  the  HHG 
radiation  is  more  intense  and  have  higher  degree  of  circular 
polarization  then  the  right-rotating  part  of  the  spectrum. 

It  is  instructive  to  analyze  the  time  domain  field  that  can 
correspond  to  the  measured  spectrum  (figure  4(a))  and  the 
degree  of  circular  polarization  (figure  4(b)).  Notably,  the 
spectral  phase  can  be  tailored  downstream  by  using  thin  metal 
films  or  multi-layer  mirrors  [70].  It  is  therefore  informative  to 
calculate  the  transform-limited  field  that  corresponds  to  the 
measured  spectrum  in  figure  4(a).  We  calculated  the  time 
domain  field  assuming  a  flat  spectral  phase  and  associate  left- 
rotating  orders  with  100%  circular  polarization  and  40% 
circular  polarization  to  the  right-rotating  HHG  orders  (i.e. 
ellipticity  of  0.21).  The  resulting  normalized  temporal  inten¬ 
sity  is  shown  in  figure  4(c)  (red  solid  curve).  It  consists  of  an 
attosecond  pulse  train  with  pulse  duration  of  ~120  attose- 
conds.  All  the  pulses  in  the  train  have  the  same  helicity. 

Importantly,  adjacent  attosecond  pulses  have  different 
degree  of  circular  polarization.  Of  course,  highly  chiral  pulses 
(high  level  of  circular  polarization)  contribute  more  to  chir¬ 
ality-sensitive  measurements  than  attosecond  pulses  with  low 
chirality.  The  dashed  blue  line  in  figure  4(c)  follows  the 
instantaneous  Stokes  parameter  S3,  which  accounts  for  the 


temporal  contribution  of  the  HHG  field  to  a  chiral  asymmetry 
measurement;  field  intensity  multiplied  by  the  degree  of  cir¬ 
cular  polarization.  To  plot  the  blue  curve,  we  project  the 
electromagnetic  field  on  left-  and  right-rotating  circular 
polarization  and  write  S3  as  the  normalized  difference 
between  the  intensity  of  the  left-rotating  field,  /Left,  and  the 
right  rotating  field,  /Right-  By  plotting  the  expected  chiral 
sensitivity,  we  note  a  new  and  interesting  feature  in  attose¬ 
cond  science.  Even  though  the  radiation  comprises  a  train  of 
attosecond  pulses,  the  contribution  to  the  chiral  asymmetry 
may  dominantly  originate  from  a  single  attosecond  pulse — 
the  most  chiral  pulse  (to  be  exact,  the  one  with  the  highest 
Stokes  parameter  S3).  In  addition,  in  figure  4(c)  we  actually 
ignored  contributions  of  the  q  =  3m  harmonic  orders  for  this 
time  domain  analysis.  In  practice,  the  suppressed  peaks  are 
left-rotating  with  30%-50%  degree  of  circular  polarization 
(see  figure  4(b)  in  [47]),  and  therefore,  may  enhance  the 
chirality  sensitivity  of  the  dominantly  chiral  attosecond  pulse 
even  further,  with  respect  to  adjacent  attosecond  pulses.  Such 
an  effectively  chirality-sensitive  isolated  attosecond  pulse 
may  lead  to  the  birth  of  chirally  selective  attosecond  science, 
where  the  dynamics  of  spin,  magnetism,  and  angular 
momentum  is  accessible  with  unprecedented  temporal  reso¬ 
lution.  Finally,  figure  4(d)  presents  a  measured  HHG  spec¬ 
trum  (recorded  in  the  same  experiment  under  similar 
conditions  as  in  figure  4(a))  with  an  even  greater  ratio  (10  to 
1)  between  left-  and  right-rotating  HHG  orders.  We  did  not 
measure  the  degree  of  polarization  for  this  spectrum  and 
therefore  did  not  analyze  its  corresponding  time-domain  field. 

8.  Conclusions 

In  this  work,  we  explored  phase  matching  of  circularly 
polarized  HHG  pumped  by  a  circularly  polarized  counter¬ 
rotating  bi-chromatic  field.  Using  two  independent  analytical 
approaches,  we  generalized  the  phase  matching  equation  to 
apply  for  a  long  bi-chromatic  pump  field  with  any  two 
spectral  components,  u\  and  u2.  In  the  first  approach,  we 
imposed  conservation  laws  for  the  energy,  momentum  and 
angular  momentum  of  the  photons  involved  in  the  HHG 
process.  In  the  second  approach,  we  analyzed  the  explicit 
mathematical  propagator  of  a  vectorial  bi-chromatic  pump 
field  propagating  in  an  isotropic  chromatically  dispersive 
medium.  We  derived  a  new  propagation  equation  corresp¬ 
onding  to  phase  matching  for  vectorial  fields.  This  vectorial 
propagation  condition  is  inherently  different  from  phase 
matching  of  linearly  polarized  HHG,  where  the  pump  and 
HHG  fields  have  the  same  polarization,  so  phase  matching 
considers  only  the  scalar  amplitudes  and  phases.  We  show 
that  the  propagation  equation  derived  from  the  propagator  of  a 
bi-chromatic  circularly  polarized  pump  reveals  a  unique 
symmetry.  It  links  three  conservation  laws  for  the  photons 
(energy,  linear  momentum  and  angular  momentum):  one  can 
derive  the  conservation  of  energy,  linear  momentum  or 
angular  momentum  from  the  phase  matching  equation 
(derived  by  the  vectorial  propagation  condition)  and  the  two 
other  conservation  laws.  We  discussed  the  symmetry  of  a 
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system  in  which  HHG  is  pumped  by  a  long  pulse  and  a 
counter-rotating  short  pulse.  By  taking  into  account  the 
effects  of  a  temporal  envelope  on  such  a  pulse,  we  found  two 
additional  regimes  for  phase  matching  and  quasi-phase 
matching,  and  showed  that  these  regimes  are  also  helicity 
selective.  In  the  last  section,  we  showed  that  the  highly  chiral 
spectrum  that  we  observed  experimentally  in  [47]  can  support 
chiral  attosecond  pulses.  We  expect  chiral  HHG  [71,  72]  and 
attosecond  pulses  to  open  numerous  opportunities  to  observe 
the  dynamics  of  chiral  phenomena,  like  magnetism  and 
diverse  charge,  spin  and  orbital  orderings,  using  all-optical 
tools  with  unprecedented  temporal  resolution  and  in  spectral 
regions  unavailable  to  circular  polarization  thus  far. 
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